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Abstract 

We propose a resolution of the renormalization group flow for the disordered Dirac fermion 
theories describing the quantum Hall transition (QHT) and spin Quantum Hall transition (SQHT), 
which previously revealed no perturbative fixed points at 1-loop and higher. The approach involves 
carrying out the flow in 2 stages, the flrst stage utilizing a new form of super spin-charge separation 
to flow to g/(l|l)Ar and osp{2\2)-2N supercurrent algebra theories, where N is the number of copies. 
This flxed point breaks the copy symmetry. In the second stage, additional forms of disorder are 
incorporated as dimension zero logarithmic operators, and the resulting actions have explicit forms 
in terms of two scalar flelds and a symplectic fermion. Multi-fractal exponents are computed 
with the result g(l — g)/4 and q{l — q)/8 for the QHT and SQHT respectively, in agreement with 
numerical estimates. 
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I. INTRODUCTION 



Disordered Dirac fermions in 2 + 1 dimensions have many important applications in con- 
densed matter physics. They are theoretically interesting since they can represent new uni- 
versality classes of Anderson localization/delocalization transitions. Perhaps the most im- 
portant is the Chalker-Coddington network model for the quantum Hall transition (QHT)[l|, 



which can be mapped onto disordered Dirac fermions 3, 



tions includes to dirty superconductors 



lattices 



[QBB 



A partial list of other applica- 



7j, and studies of hopping models on bipartite 
More recent applications are to graphene ^, [l^ , where the Dirac fermions are 
present from the start. The possible universality classes of disordered Dirac fermions were 
classified according to their discrete symmetries in ll|. The latter classification contains 13 
classes and is thus a minor refinement of Altland-Zirnbauer's classification which does not 
assume the Dirac structure fl2|. 

A number of new theoretical techniques have been developed over the last decade to 



13 



14 



15 



3, HQ, 



19|. For the 
acking. 



study these problems; a partial list includes for instance 

most part, a proper understanding of the critical points for generic disorder is still 
A notable exception is the spin quantum Hall transition (SQHT). Its network model 
also be mapped onto disordered Dirac fermions 20|] . Remarkably, the equivalent spin chain 



2l[ | can 



was mapped onto 2D classical percolation by Gruzberg, Ludwig and Read[22|, |23|, |2J], and 
this leads to the exact knowledge of the correlation length exponent I'perc = 4/3 and density 
of states exponent p{E) ~ E^^'^. 

For the QHT, one should also mention the replica sigma model approach of Pruisken 25]. 
Although it appears to have the right ingredients as outlined in j26|, it has proved too difficult 
to solve thus far, so it remains unknown whether it really does have the correct critical 
point. There is also the later proposal of Zirnbauer 27| which uses supersymmetry. Based 
on symmetry and various other requirements the critical point for the QHT was proposed 



to be described by a sigma mode 
model was further studied in 



of WZNW type based on the supergroup PSL{2\2). The 
. The main problem with this proposal is that the level 
k of the PSL{2\2) WZNW model is an exactly marginal perturbation so that the model 
actually has a line of fixed points depending on k. This would lead to the prediction of 
non-universality in the QHT, which is contrary to the numerical evidence. (For a recent 



review, see 



29 



32l|.) It was pointed out recently by Tsvelik that the value k = 8 gives 
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very reasonable exponents [33|. (The identical exponents were actually already speculated 



m 
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371].) Unfortunately, it was noted that there are no known constructive arguments 



leading to A; = 8 based on the PSL{2\2) approach. In the work we present here, PSL{2\2) 
will not play a role, but rather the simpler superalgebra gl{l\l) will be central, and we will 
describe a precise mechanism for obtaining higher integer levels k based on a new form of 
super spin-charge separation [35|. 

Due to an extensive effort over the last few decades, vast classes of conformally invariant 
2D critical points can be constructed [s^ and many of these mathematical constructions can 
be extended to theories with supergroup symmetries. However in the study of disordered 
Dirac fermions, it is important that the possible critical point is supposed to be reached by 
renormalization group flow in the effective disorder- averaged effective field theory. There 
are a very limited number of known mechanisms for obtaining a fixed point from a renor- 
malization group flow, and unless disordered systems depend on some new mechanism, it is 
helpful to identify the known mechanisms guide: 

(i) Non-linear sigma models like the 0{N) model. Due to the Mermin- Wagner theorem, 
these models are only critical for —2 < N < 2. 

(ii) Non-linear sigma models with topological terms, the primary example being the 0(3) 
non- linear sigma model with 6 = tt. 

(iii) Certain relevant perturbations that induce flows between minimal models. 

(iv) Decoupling of massive degrees of freedom using spin-charge separation in the case of 
marginal current- current interactions. Here the primary example is the Id Hubbard model 
which has SU{2) (g) SU{2) symmetry. One SU{2) sector is marginally relevant, the other 
marginally irrelevant. Thus one SU{2) is gapped out in the flow, and the fixed point is the 
SU{2) current algebra theory (WZNW model) |31j. 

For the sake of comparison with our work, Pruisken's model is based on the scheme (ii), 
whereas Zirnbauer's model is mathematically constructed directly at the fixed point so it is 
unknown under which renormalization group flow scheme it can be realized. In the approach 
pursued in the present work, disorder averaging of Dirac fermions is known to yield marginal 
current-current interactions so naturally the mechanism for obtaining a fixed point will be 
(iv), i.e. based on super version of spin-charge separation. 

The starting point of the present work is the detailed disordered Dirac fermion theories 
for the Chalker-Coddington network model for the QHT and its variant for the SQHT. 
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Based on extensive numerical evidence 



321 ]. there is no doubt that these specific models 



have a critical point. Performing the disorder averaging using Efetov's supersymmetric 
method 36|] is known to lead to marginal anisotropic left-right current-current interactions 



of the underlying super current algebra. The problem with such marginal perturbations is 
that they typically do not have perturbative fixed points at finite values of the coupling 
constants. In particular, the coupled beta functions do not have any non-trivial zeros at 
one loop. Higher loop corrections to the renormalization group (RG) beta functions were 



computed for the network models in SJ, l37| based on the general proposal in 



and also 



did not reveal any perturbative fixed points. These analyses were nevertheless useful for 
understanding whether any new couplings were generated under RG. In these studies it is 
significant that the couplings fiow to a singular point in a finite RG time, which suggests an 
incomplete resolution of the fiow rather than the lack of a fixed point. It was pointed out 
that the higher order beta functions are possibly not exact due to some contributions that 
were missed at 4-loops 39|| , however it seems unlikely that this could resolve the issue in a 
constructive manner. 

Since a perturbative fixed point of the beta functions is unlikely, one must identify the 
correct non-perturbative mechanism that singles out the expected fixed point. In this paper 
we propose to resolve the RG fiow in two stages. We first focus on the important symmetries 
of the N-copy theory before disorder averaging, i.e. we identify the relevant symmetries that 
are present for any realization of the disorder. This leads to the special role of the symmetries 
corresponding to the current algebras gl{l\l)N ® su{N)o for the QHT and osp(2|2)_2Ar ® 
sp{2N)q for the SQHT. (Our nomenclature is that refers to the current algebra for the 
finite (super) Lie algebra g at level k. ) The disorder averaged effective actions have several 
couplings which correspond to the strengths (variances) of the various random potentials. 
Rather than study the simultaneous flow of all couplings, in the flrst stage we set some of 
the couplings to zero and carry out the RG flow for a subset of the couplings corresponding 
to the above symmetries. A new form of the super spin-charge separation obtained in 35|] is 
then used to argue that in the flrst stage one flows to the flxed point gl{l\l)N for the QHT. 
This result indicates that the locality constraints studied in 13|] for the gl{l\l)k theory, 
which led to k being an inverse integer, are too restrictive. For the N-copy SQHT, the 
an alog ous flow is to to osp{2\2)_2N- This kind of flow for 1-copy of the SQHT was studied 
in 171] where it was viewed as a flne-tuning of the model. In this paper our point of view 
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is that the first stage of the RG flow identifies the proper degrees of freedom that are the 
most important for the actual critical point. 

In the second stage of the flow, we restore the additional kinds of disorder that were 
possibly initially present in the model as additional relevant perturbations. The possible 
operators which appear in this second stage are dictated by the quantum numbers of the 
original fields and the super spin charge separation. Another new aspect of the present work 



is that we use the results in 



35l | to explicitly construct the operators corresponding to the 



additional kinds of disorder. In particular, (yf/(l|l)fc at any level k has a simple free field 
representation in terms of two scalar fields and a symplectic fermion. The additional kinds 
of disorder correspond to logarithmic operators of scaling dimension zero. For the QHT 
one obtains a gl{l\l) generalization of the sine-Gordon theory, where A^, the number of 
copies, appears as a coupling. It can also be viewed as radius of compactification R = ^/N. 
The important feature of these kinds of perturbations is that they do not drive the theory 
to a new fixed point, but rather just lead to logarithmic corrections to the correlation 



functions 



35 



401], and this explains why for example the osp(2|2)_2 current algebra contains 



the correct exponents for the SQHT. 

As models of disordered Dirac fermions, the QHT and SQHT are not so different in their 
formulations, and if the methods are general enough, they should be subject to the same 
kind of analysis. It is therefore very instructive to work out both cases in parallel, since 
some exact results are known for the SQHT. This also avoids idiosyncratic proposals for 
special cases. All the remaining sections of this paper have subsections treating the QHT 
and SQHT cases. 

Our results are presented as follows. In sections II and HI we review the definitions of 
the models and the supersymmetric method for disorder averaging, introducing a convenient 
notation to deal with the profusion of fields in the N-copy theories. The symmetries of the 
models for any realization of disorder are studied in section IV. In section V we consider a 
subgroup of these symmetries that commutes with the permutation of the A^ copies, which 
leads to 5'/(l|l)Ar and osp{2\2)^2N- The properties of these super current algebras that we 
need are reviewed in section VI. In section VII we describe our 2-stage strategy for resolving 
the RG flow. Since we focus on symmetries that are present for any realization of disorder, 
the analysis does not depend strongly on any assumed distributions of the random poten- 
tials. On the other hand the arguments rely strongly on super spin-charge separation and 
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some simple 1-loop beta function arguments concerning the marginal relevance/irrelevance 
of operators in the disorder averaged effective action. The additional perturbations in the 
second stage of the RG flow are constrained by the quantum numbers of the fields after 
gapping out the su{N)q and sp{2N)o "copy" symmetries. Typically one obtains pertur- 



This results in 



hamiltonian[2, l3| 



bations by logarithmic operators which were explicitly constructed in 
some relatively simple lagrangians involving the two scalar fields and a symplectic fermion. 
Under some assumptions, the multi-fractal exponents are computed in section VIII. Our 
results agree favorably (within about 1%) with the numerical simulations in 4ll. |42|. l43l. |44|. 
In section IX we discuss the localization length exponents. 

II. DEFINITION OF THE MODELS. 

A. Chalker-Coddington network model. 

The Chalker-Coddington network model can be mapped to the following 2d 

V + M -td, + AA 
H={ 1 
-id, + A, V-M J 

where z,z are euclidean light-cone coordinates, z = {x + iy)/^/2,z = z* with x,y the 2d 
spacial coordinates. A^^z is a u{l) gauge field, A^ = (A^ + iAy)/^/2. The hamiltonian is 
hermitian if A^, Ay, V, and M are real. The hamiltonian is first order in derivatives and 
operates on a 2-component wave-function. It thus corresponds to a universality class of 
disordered Dirac fermions, class A=GUE in p^], or class according to the more specific 
classification i n All the potentials A,M,V depend on x,y and are random variables. 
The model in 16| on the other hand is in the chiral GUE, i.e. class 2. 

B. Spin network model. 

The network model for the SQHT is also a model of diordered Dirac fermions, but in 
class C [l^ (class 4_ according to[ll|). The hamiltonian is 

'2a-a + M -idz + Az\ 
H={ (2) 

-id. + A, 2a- a- M 
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where a are Pauli matrices and A is an su{2) gauge field, A = A - a. The hamiltonian thus 
operates on a 4-component wave-function. Again, all the potentials a, M, and A depend on 
X, y and are random. 



III. SUPERSYMMETRIC DISORDER AVERAGING. 



Since the hamiltonians describe non-interacting fermions, the disorder averaged correla- 
tion functions can be studied with Efetov's supersymmetric methods^. 



A. QHT 

Let us denote the 2-component wave-functions as follows: 

^=(|^^j, r = {^-,i^-) (3) 

For simplicity, let us refer to all the disordered potentials simply as "V . The Green functions 
can be defined with respect to a functional integral with the action 

S{i^;V)=i I ^rH{V)^P (4) 

Above, the fields are taken to be fermionic. 

The supersymmetric method is a trick to cancel the fermionic determinant Z{V) where 
Z is the partition function at fixed disorder V. One introduces ghost partners /3, to the 
ip^s and considers the action 

Ssusj = Si^; V) + SiP; V) (5) 

where S'(/3; V) is identical to S{ip; V) but with the replacement ip ^ [3. The /3-fields are 
bosonic. The effective action upon disorder averaging is then defined as 

Q-SMP) = j DVV[V]e-^'^'^ (6) 

where 'P\y] is the probability distribution of the random potentials. If 'P\y] is taken to 
be gaussian, then S'efr contains quartic interactions among the fermions and ghosts. As we 
will see, many of our arguments are independent of the specific form of these probability 
distributions. 
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In order to clearly display the symmetries of (Sgusy and Ses, it will prove convenient to 
introduce the following notation. Let ^± denote 2-component fields built of out ■0±, /3±, and 
similarly for ^-|-: 

*± = (^±,/3±), *± = (?±J±) (7) 
The index that runs over the two components of ^± will be denoted as "r", r = 1,2: 

We will also be interested in computing disorder averaged moments of correlation func- 
tions. To compute averages of iV-th moments such as 

we need to introduce A^-copies of the models. Namely, we introduce fields a — 1,..,N, 
so that the complete set of fields is and \E'2", r — 1,2. Thus refers to 2A'" different 
fields. 

At a fixed realization of disorder, 5'susy can be expressed in the compact form: 

- iM(x) (^_*+ - *_^+) ] (9) 
where for example "^^"^^ = Ylra^-"^+"- 



B. SQHT 

For the spin-network model, one needs to introduce an additional su{2) index "i" and 
consider fields i.e. there are 4A'" fields in ^+ for example. The action is then 

-^susy = / ^[ ^-(9. - iA,{x))^+ + - - ia{x) ■ + 

-iM{x) {^-^+-^-^+) ] (10) 

Above, the Pauli matrices, including the ones in A, operate on the index i so that for example 
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IV. SYMMETRIES AT FIXED DISORDER. 



A. QHT 

First consider all disordered potentials set to zero in 5'susy. The result is a free conformal 
field theory of Dirac fermions and ghosts which has total central charge equal to zero: 

72. ^ 



^free = / ^ E (^-^^^+ + ^-^^l + /^-^^/^+ + f^-^f^l) (H) 

The two point functions are 

(z^_(^)z^+H) = (z^+(z)V^_H) = {P^{z)p.{w)) = -{(34z)f3^{w)) = (12) 

z — w 

and similarly for the right-movers, {ip _^{w)) = l/{z — w), etc. (For a review of 2D 
conformal field theory see 30|, |45|].) In the sequel we will not display the right-moving 
counterparts if they are the obvious duplications of the left. 

For notational simplicity let us group the r, a indices of into a single index a and 
refer to these fields as a = 1, ..,2N. The extra minus sign in the above two point 

functions can be accounted for by introducing a grade [a] = for bosonic components and 
[a] = 1 for fermionic ones. One then has 

^^(z)^^(O) ~ - ^^(2)^^(0) ~ - (_)H+i5»'' (13) 

The complete set of chiral currents are then 

These currents generate an osp{2N\2N)k super-current algebra at level k = 1, and this 
represents the maximal symmetry without disorder. Our conventions for all the current 
algebras that appear in this paper are presented in Appendix A. 

Without disorder the symmetry is actually the sum of left and right, osp{2N\2N)i © 
osp{2N\2N)^ since the theory is conformal. In the presence of disorder the conformal 
symmetry is broken and one does not have the full current-algebra symmetry. However one 
can study the global left-right diagonal symmetries generated by the charges 

Q = / ^ A^) + / ^ m (15) 



where J, J are the left / right moving currents. These conserved charges are always associated 
to a set of left-moving currents for a current algebra and in the sequel this correspondence 
is implicit. 

With disorder, the maximal osp(2A^|2A^) symmetry is broken to something smaller. Con- 
sider the transformation ijj ^ ijj + 6ijj which acts left-right diagonally: 5?/'" = (3"' , 6(3^ = 
— = l3_ , 6(3_^_ = and is zero on all other fields, where a, a' are fixed copy 

indices. All of the operators \E'_\1'+, \1'_\E'+, ^_\1'+ and are invariant under this trans- 

formation. There is another symmetry of this type with + ^ — . The left-moving currents 
that generate these two symmetries are 5"^'° = . The charges for these fermionic 

symmetries are nilpotent, = 0, and the symmetry they generate can thus be thought of 
as a BRST symmetry. Namely, the disorder dependent part of S'susy can be written as 6X 



for some X, and its invariance is a consequence of 5^ = 14| . 

Consider other nilpotent symmetries with 5^p" = i.e. that flip the u{l) charges. 
These correspond to the currents S^'°' = ip^f3^ ■ One finds in this case that due to fermionic 
exchange signs, the operators \E'_\I^_|_ and ^'.^'.i- are not invariant. The only invariant is the 
combination (l'_\l>+ - \1>_"5'+). 

Thus, examining the action eq. iQ, one sees that the only nilpotent symmetries are 
those corresponding to the diagonal left-right symmetry which corresponds to the left-moving 
currents S*^'" . The operator product expansion (OPE) of these currents closes on the super- 
current algebra gl{N\N)k=i- (See Appendix A.) In other words, at any fixed realization of 
the disorder, the model has a global gl{N\N) symmetry corresponding to the current algebra 
gl{N\N)i. We will refer to this symmetry as the BRST symmetry. 



B. SQHT 

Since there are twice as many fields in the SQHT, the maximal current algebra symmetry 
with zero disorder is osp{AN\4N)i. Repeating the analysis above for the QHT, one sees that 
there are nilpotent symmetries generated by the currents S"^"' = ± X]i=i 2 V'i.j/^^'r 

A basic result we will use repeatedly is the following. Given two copies of the same 
current algebra with currents at level ki and J2 at level k2 which furthermore commute, 
[Ji{z), J2, (w)] = 0. Then J° = + J2 satisfies the current algebra at level ki + k2. Thus, 
since the su{2) indices i are summed over in S±, these currents close on gl{N\N)k=2 since 
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each copy has level 1 and the levels add. 

The gl{N\N) symmetry is actually enlarged due to an additional nilpotent symmetry. 

Introduce the matrix e which acts on the su{2) indices i: e = ^ ^ ^ Using ea = —a^e 

and = —1, one can verify that S'susy is invariant under the left-right diagonal symmetry 
corresponding to the left-moving currents S*^"' = ijj^ ^P±- It is important to note that this 
symmetry would not be valid if there were additional type of disorder, or if the gauge 
field A contained a U{1) component. For = 1 copy, the currents S±, S± close on the 
osp{2\2)k current algebra at level k = —2^^. For A^-copies this symmetry is promoted to 
the BRST symmetry osp(2A^|2A^)_2. 



V. PERMUTATION INVARIANT BRST SYMMETRIES. 

The BRST symmetries discussed in the last section are rather large since their dimension 
depends on the number of copies N. Furthermore, the current algebras are only moderately 
interesting as possible critical points; for example the gl{N\N)i theory has only integer 
scaling dimensions at level 1. In this section we constrain the possible fixed point further 
by considering permutations in the number of copies. 

Let Vn denote the discrete permutation group for elements. The actions S'susy possess 
this symmetry where N is the number of copies. It is natural then to make the hypothesis 
that a possible fixed point also has the permutation symmetry. The BRST symmetries of 
the last section do not commute with Vn, however there is a sub- algebra that does, which we 
will refer to as the Vn invariant BRST symmetry. In section VII we will provide arguments 
based on super spin-charge separation that indicate how a fixed point with this restricted 
symmetry can arise under RG flow. 

A. QHT 

For the QHT the generators that commute with Vn are 

H = J2nr-, J = 5^/3^/3^, S^ = ±J2^Zf3^ (16) 

a a a. 

The above currents satisfy the gl{l\l)k current algebra at level k = N\ 

H{z)H{0) ~ A j(^)j(o) ^ -A 
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Hiz)S±iO) ~ J(^)S±(0) ~ ±- S± (17) 

Z'^ z 

It will be important to determine any additional continuous symmetries that commute 
with the Pjv-invariant BRST symmetry gl{\\\)jq. There is obviously an su{N) symmetry 
which mixes the copies. Let L'^,Lp denote the su{N) currents in the separate sectors and 
L"- their sum: 

^; = CC'<, l; = /3%,,/3:', l'^ = l; + l; (is) 

where is the N x N dimensional matrix representation of the vector of su{N). The 
currents satisfy su{N)i, whereas the Lf^ satisfy su{N)^i. Therefore the total currents 
satisfy su{N)q at level k = 0. In summary, the symmetries that will play a significant 
role in the sequel is gl{l\l)iy © su{N)q and these two current algebras commute. 



B. SQHT 

For the SQHT the Vn invariant BRST symmetries correspond to the currents 

H = 13+13. , J = V- , J±=^P±e^P± (19) 
S± = i^±f3:f, S± = 'ilj±e(3± 

where ip+ip- = a'^'+,i'^-,i 4'±^'^± = Si j a "^±,4 ^ij "^ij; 6tc. The above currents 
satisfy osp(2|2)_2Ar. 

The S± and J± currents are invariant under iIj± — ^ MiIj±, (3± M(3± where M is a 2N 
dimensional matrix satisfying M*(e ® 1)M = e®l. M is thus an element of Sp{2N). The 
currents satisfy sp{2N)q. Since S± and J± close on osp(2|2)_2Ar, the sp{2N)q commutes with 
osp{2\2)_2N 

For = 1, sp{2) = su{2), and this su{2) is the original 5^(2) symmetry of the spin 
network model. We wish to emphasize that here the sp{2N)q symmetry is a property of the 
A^-copy theory, which is to be contrasted with different models that have a random sp{2N) 
gauge field from the very beginning in the 1-copy theory, for example the models inj^. 
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VI. THE gl{l\l)k AND osp{2\2)k SUPER CURRENT ALGEBRAS. 



9l{ 
on 



In this section we summarize the main results we will need for the current algebras 
l) h an d osp{2\2)k. For gl{l\l)f: we mainly summarize our recent work 35|], which builds 



47, 



48| . The osp{2\2) results are based on the workjlSj. 



A. gl{l\l)k 

We will need the Sugawara stress tensor T{z). The algebra gl{l\l) has two independent 
quadratic casimirs: 

C^ = J^-H^ + S+S- - S-S+, C'^ = {J- Hf (20) 

where it is implicit that the above operators are the zero modes of the currents. The stress 
tensor is fixed by the condition T{z)J"'{Q) ~ J"(0)/2;^, which requires it to be built out of 



both casimirs 



49|: 



nz) = {J' -H' + S^S^ - S^S^) + ^(J - HY (21) 

For any level k there exists a free field representation in terms of 2 scalar fields and a 
symplectic fermion. The free fields have the action 

2 

S = ^ I J2 {Vabd.rd,^' + eatd.x'^d.x'') (22) 

^ ^ a,b=l 

where 

and (9^ = 2dzdz- The x fields are Grassman: (x*^)^ = and have Virasoro central charge 
c = —2, so that the total central charge is zero. Note that the metric for the bosonic fields 
has indefinite signature, which will turn out to be important. The equations of motion imply 
that the fields can be decomposed into left and right moving parts: 

cf>''iz,z) = (24) 

X^iz^z) = x'{z)+riz) 

Above, the bold face signifies local fields. The two point functions are 

{r{z)cl>\w)) = -v^Hog{z-w), {x''{z)x\w)) = -e'^Hog{z-w) (25) 
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(Our conventions are t]"''^ = rjab, ^""^ = ^ab-) Exponentials of the bosons have the conformal 
dimension: 

^ (^g^(a0l+602)j ^ (26) 



It is straightforward to verify the following representation of the OPE's in eq. ([Tj 

H = iVk,d,(f)\ J = tVkd,(f)'^ (27) 

In the sequel, where there is no cause for confusion, we will simply write d(j) for dz<f){z). 
The algebra gl{l\l) has 2-dimensional representations where H = diag(/i, /i — 1) and 



J = diag(j,j — 1) which will be denoted as {h,j). (We follow the conventions in |35|.) 
These are so-called typical representations when h ^ j. Primary fields associated with these 
representations have conformal dimension 

^i^.) = + (28) 

The basic fields 4'±, P± are in the fundamental representations: 

(^+,/3+)^(l,0), (/?„,^_)^(0,1) (29) 

and have scaling dimension 1/2 when k = 1. 

The tensor product of two typical representations is 

{hi,ji)®{h2,j2) = {hi + h2,ji+j2)®{hi + h2 - l,ji +j2 - 1), {hi + h2 ^ ji+32) (30) 

When hi + h2 = ji + 32 the tensor product gives a new reducible but indecomposable 
4-dimensional representation denoted {h)(4^y. 

(hiJi) ® (/i2, J2) = {hi + h2- 1)(4) (31) 

These representations have A(/j)j^j = 0, however they are logarithmic since the casimir C2 is 
not diagonal. 

The vertex operators V(^h,j) can be explicitly constructed in the free field theory, and 



require the twist fields in the symplectic fermion sector 35|]. As for the spin fields of the 



Ising model, the twist fields modify the boundary conditions of the fundamental field x'- 

X^(e2-^^)/XA(0) = e-2-^\i(^)^A(0) (32) 
X^(e2-z)/x,(0) = e2-\2(z)/x,(0) 
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The properties of these fields were studied in[50|. It is clear from the above equation that 
2ti\ is a phase and is restricted to — 1 < A < 1. We also need the doublet of twist fields a^, 
which arise in the OPE: 



The conformal dimensions of the twist fields are 



A(A-l) 



A 



(x) 

A ' 



A 



(x) 



Afa?l 



a{x) 



(33) 



(34) 



The vertex operator for V(^h,j) requires twist fields with X = {h — j)/k. For h > j: 



whereas for h < j: 



{h-j) 



1/4 



-/iA e'(^'^'-^'^')/^ 

^2 ^i{{h~l)<^>^~(j-l)4>^')/^/k 
A 



V, 



{h,j) 



1 ^i{h4>'^~j<l>^)/Vk 



A 



h-j 
k 



h-j 



(35) 



k 



(36) 



U - hfi' 

The vertex operator V/^h)^^^ for the representation (/i)(4) is constructed from the logarithmic 
field eabX'^x''- 



(4) 



/ ^1 gi(h+l){0l-</,2)/yfc , 
^ gi/l(0l-02)/v/fc 
J_^1^2 pih{4,^-,j>^)/Vk 



(37) 



y ^2 gi(/i-i){</,i_</,2)/v^ y 

The two middle fields form a logarithmic pair with A = 0. It is important that the above log- 
arithmic field has a simple and explicit construction in the second-order symplectic fermion 
theory; this is not transparent in the minimal model description of c = —2, nor in the 
first-order description. For a review of logarithmic conformal field theory, see 

A closed operator algebra is obtained when k is an integer and the spectrum of fields 
V(/ij) is restricted to h,j integers satisfying 



k < h — j < k 



(38) 



This operator algebra can be viewed as generated by OPE's of the fundamental vertex 
operators V(i_o) and V(o,i). For instance: 



V(i,o>(2;)V(o,i>(0) ~ -17^ 



(39) 
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B. osp{2\2)k. 



The finite dimensional representations of osp(2|2) can be labeled by the su{2) with gen- 
erators J, J± and by the u{l) charge H. The typical, irreducible representations will be 
denoted as [6, s]"^"^ where s G {0, |, 1, |, ....} is an su{2) spin and h = H/2. These represen- 
tations are 8s dimensional. In order to describe their 5^(2) ®u{l) decomposition, let [b, sY^ 
denote the 2s + 1 dimensional representation with J/2 = = — s, — s + 1, s and H = 2b. 
The generic decomposition is 

[b, sY'P = [b, s]^" © [6 + i, s - i]^" © [6 - i, s - i]^" © [6, s - 1]^" (40) 

As for 5f/(l|l), there are atypical, indecomposable but reducible representations at 6^ = s^. 
The simplest is 8 dimensional and arises in the following tensor product 

[0, lY'P © [0, '^Y'P = [0, 1]°^^ © [8]°^^ (41) 

The [8]°'P has the quantum numbers of [i, © [-|, |]°"^. 
The stress tensor is built from the single quadratic casimir: 



T ^ 

-'osp(2|2) 



2(2 - k) 



(42) 

and the typical representations with 6^ ^ have conformal dimension 

2(s^ - b^) 

^lb,s] - 2-k ^^^^ 
At A; = —2 there is a free field representation with the same field content as in eq. (!22|) \^ . 



This can be derived from the gl{l\l)2 embedding 35|]. In fact: 



Tosp{2\2)^2 = Tgl(l\l)2 = Tgl(l\l)_2 (44) 



VII. DISORDERED CRITICAL POINTS: SUPER SPIN-CHARGE SEPARATION 
AND THE RENORMALIZATION GROUP. 

A. General strategy. 

If the random potentials are taken to be gaussian distributed, e.g. VlV] = 
exp(— ^^j^ J (Px V{xY), then the functional integrals over all the random potentials can 
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0^ = Y.d^,.r.f (46) 



be performed and one obtains the general form 

Scff = 5free + / ^ E 3aO\x) (45) 

where Si^ec is the free action of the iIj±^ [3± fields eq. ffTTj) . (/yi are variances that measure the 
strength of the disorder, and are marginal operators. The operators can always be 
expressed as left-right current-current perturbations, i.e. they are of the form 

a,b 

for some bilinears c?^^, and the currents J" are those of the maximal symmetry, osp{2N\2N)i 
for the QHT and osp(4A^|4A^)i for the SQHT, and are bilinears in the fields il'±y Since the 
models have the BRST symmetries for any realization of disorder, the operators must be 
BRST invariant, i.e. Ses has a gl{N\N) symmetry in the QHT and an osp{2N\2N) symmetry 
for the SQHT. The perturbations can thus be viewed as anisotropic interactions of the 
maximal current algebra that are BRST invariant. 

As discussed in the introduction, the perturbative RG for the simultaneous flow of all the 
couplings qa does not reveal a fixed point. In order to resolve this difficulty, we propose to 
perform the RG in two stages, with special attention paid to the symmetries that exist at 
any realization of disorder. 

We will need the following general property. Consider two commuting current algebras 
Qa and Qb with currents Ja, Jb- Furthermore, let us suppose that the stress tensor for a 
given conformal theory separates as follows: 



-eft 



Tg^ + Tg, (47) 



Consider the perturbation of the conformal field theory by left-right current- current pertur- 
bations: 

r (fx — — 

^ = 5cft + / — {gA Ja-Ja + 93 Jb ■ Jb) (48) 

where J ■ J is the invariant built on the quadratic casimir. Since the currents commute, the 
RG beta-functions decouple; to 1-loop the result is: 

^ = Cf9l ^ = Cfgl (49) 

where i is the logarithm of the length scale and is the casimir for the adjoint rep- 
resentation of the finite dimensional part of Qa- Let us suppose that the physical regime 
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corresponds to positive gA,B- If is positive, then the couphng Qb is marginally rele- 
vant and the flow is to inflnity. This is a massive sector and the Qb degrees of freedom 
are "gapped out" in the RG flow. If is negative, then the coupling is marginally 
irrelevant. This results in the flxed point characterized by the Qa current algebra. If the 
original conformal fleld theory corresponds to the current algebra ^max, then the flxed point 
may be viewed as the coset Gmax/GB- This scenario was proposed for generic flxed points of 
marginal current-current perturbations in 3J], however what was missing in the argument 
was the spin-char ge s eparation. The resulting coset is a somewhat trivial example of the 
GKO construction 53|] because of the separation of the stress tensor. 

In the flrst stage of the RG flow we consider only the disorder couplings for the Vn 
invariant BRST symmetries and the additional "copy" symmetries. The essential ingredient 



here is the spin-charge separation 17|, 



35 



46] . In the flrst stage we will obtain a critical point 



corresponding to the Vn invariant BRST symmetry. In other words, in the flrst stage of the 
RG flow we identify the massless degrees of freedom that are most important at the critical 
point. In the second stage of the RG flow, we reintroduce the other kinds of disorder as 
further relevant perturbations of these massless degrees of freedom. 



B. QHT 

The Vn invariant BRST symmetry is gl{l\l)i\i and the copy symmetry that commutes 



with it is su{N)q. Remarkably there exists the following spin-charge separation 



35]: 



1 ^ 

tL''"'" = -2 J2(r-d.n + ^-d^n) = T^im... + TsuiN), (50) 

a=l 

where the stress tensors Tgi(^i\i)^ and Tsu{n)o are the Sugawara stress tensors for interacting 
theories. Simple checks of the above result are as follows. First, all of the stress tensors 
have c = 0. The free theory contains AN flelds ip'^, f3'^. Under the gl{l\l) (S> su{N) they 
transform as ((1,0) © (0, 1)) (8> [vec] where [vec] is the A^-dimensional vector representation 
of su{N). The later has conformal dimension Asu{N)k = 2N{k+N) general level k, whereas 
A(i^o) = ^(o,i> = 2p'- Oiie sees that the dimensions add up properly: A(i_o) + As^i^j^^^ = 1/2, 
as is appropriate for the free ip±, [3± flelds. 

In the flrst stage of the RG flow we consider the action of the form eq. (HHj) where 
Ga = 5'^(l|l)Af and Qb = su{N)q. For su{N), C^'^^ > 0, and it is gapped out in the flow. For 
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gl current perturbations the situation is somewhat more subtle because there are two 

quadratic casimirs 



161 ] ■ Consider 

S = Stoc + [—[g {JJ-HH + S+S^ - S^S+) + g'{J - H)(J - H)] (51) 
Then the 1-loop beta function for g is zero, whereas dg'/di = —g^^^. (Here we fixed 



the sign by setting some couphngs to zero in the more general result in [3J]. ) Therefore 
these gl{l\l) current interactions are irrelevant. It is important that this is in contrast to 
the situation for the model in the Gade-Wegner universality class which has the same beta 
function up to a sign[l^, dg'/di = +g'^, which implies the disorder in that case is marginally 
relevant. This difference in sign is a consequence of the detailed form of the hamiltonian 

n 

in pj| which acts on a 4-component wavefunction rather than two. The sign of the beta 
function can also be flip ped by considering an imaginary gauge potential. The higher loop 
corrections computed in 3J] do not alter this picture. Thus in the first stage of the RG flow, 
one flows to the fixed point gl{l\l)]\[. 

It is important to point out that one feature of this scheme for ariving at the fixed point 
gl{l\l)N is that the critical exponents now depend on A^, whereas in the original gl{N\N) 
invariant theory they were independent of A^. This is due to the fact that we have gapped 
out the SU{N)q sub-algebra of the gl{N\N)i. On the other hand, the common assumption 
in the literature is that the critical exponents should be independent of A^, even though an 
RG scenario that achieves this has not yet been proposed. We will return to this issue where 
it must be faced in section VIII on the multi-fractal exponents. 

Additional kinds of disorder can now be incorporated as follows. The original currents for 
the maximal symmetry osp(2A^|2A^)i, eq. iHM can be classified according to the 5'/(l|l)Ar ® 
su{N)q. Since the RG flow in the first stage gaps out the su{N)q, what remains are the 
gl{l\l)iy degrees of freedom. The possible gl{l\l)N representations the residual currents 
fall into follows from the fact that they are bilinears in the fields ip±, f3±, and the latter 
correspond to the (1,0) and (0, 1) representations, eq. ( l29i) . Using the tensor products 

(1,0)0(1,0) = (2,0)0(1,-1) 

(0,1)0(0,1) = (0,2) ©(-1,1) (52) 
(1,0)0(0,1) = (0)(4) 



one sees that the currents transform in the (2, 0), (1, —1), (0, 2) and (—1, 1) representa- 
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tions, whereas H"^^ transform in the (0)(4). Note that (0)(4) has the same quantum numbers 
as the adjoint of gl(l\l), and represents what is left of the gl{N\N) after gapping out the 
su{N)q. The conformal dimensions of these representations are 

A A 2 + AT 

A(2,o> = = 

2- N 

A(o,2> = A(i,_i) = (53) 

and A(o)(4) = 0. 

Let ^r{z,z) denote the gl{l\l) invariant local field associated with the representation r 
of 5f/(l|l) with scaling dimension 2Ar. It can be expressed as a product of left-right vertex 
operators Vr ■ Vr- For general N, the fields $(2,0), ^(0,2) and are expressed 

in terms of the twist fields fi\ and with A = 2/N, 1 — 2/N. Explicit expressions can be 
found in [35]. The field $(o>(4) on the other hand requires only the symplectic fermion and 
bosons, and for < 2 it is the most relevant operator. If we keep only the most relevant 
operator, then we should consider 

^ = ^gimN + 9 J ^ *(0>(4) 

= / ^ ( E ^afe d.rd.^j^" + eat d.x^'d.x' + 9 x'x' cos ((0^ - 

\a,b=l 

where here, and henceforth, 0, x are the local fields ^, X- 

The case of = 2 is distinct since both $(0,2) and have A = which is degenerate 

with the dimension of ^(o)(4)ji_so that the latter is no longer the most relevant operator. The 
explicit forms at = 2 are 35l |: 




*(!,_!) - $(o,2> = ^xV COS ((01 + 0')/y2) + ^d^X%x'){xV) cos(y202) (55) 

($(0,2) only differs from $(0,2) by some fermionic exchange signs.) For N > 2 the opera- 
tors ^*(o,2) become more relevant than $(o>(4) and may need to be included as additional 
perturbations in eq. fl5^ . 

The important feature of logarithmic perturbations such as eq. fl54l) is that we believe 
they should not drive the theory to another fixed point, but rather just give logarithmic 
corrections to correlation functions. General arguments were given in 40] for marginal 
logarithmic perturbations. Here the logarithmic perturbation has dimension zero and is thus 
strongly relevant so the arguments in 40| do not necessarily apply. However arguments based 

20 



on the fact that the perturbation has exactly dimension zero were given in [35|] supporting 
the idea that one is not driven to a new fixed point. For the concrete model eq. fl54l) this 
is easy to see since, due to the indefinite signature of the scalar fields, the OPE of the 
exponentials is regular: 

This implies that in conformal perturbation theory the perturbation behaves like a mass 
term x^X^^ ^ind does not lead to contributions to the beta function at any order in pertur- 
bation theory. A contribution to the beta function for g would require a singular term in 
the OPE of the form $(o)(4j (x)$(o)(4j (0) ~ ^(o>(4)) however there is no such term. Although 
we cannot give any stronger arguments at this stage, in the sequel we will adopt the work- 
ing hypothesis that dimension zero logarithmic perturbations essentially do not change the 
critical exponents. 

C. SQHT 



The version of super spin-charge separation we need was proven in 46] : 

^fr^ ^"^^ = ^osp(2|2)_2jv + Tsp{2N)o (57) 

Here the check of the scaling dimensions goes as follows. The 8A^ fields ip'^^, (3^^ transform 
as [0, i]"'^^® [vec] under the osp{2\2) ^2N ® sp{2N)o where the vector representation of sp{2N) 
is 2N dimensional. For general k the latter has dimension Asp(^2N)k = 4{k+N+i) ■ level 
k = —2N from eq. (143|) one has A°q^i-^ = ^n+i) ■ Again one has + ^sp{N)o = 1/2, as 

required. 

Repeating the same arguments as for the QHT, since C'^p(2Ar) ^ ^ ^osp(2\2) ^ 0' 
first stage of the RG fiows takes us to the current algebra osp{2\2)^2N- For = 1 this RG 



fiow was studied in greater detail in 17|, |37[^ where it was viewed as a fine-tuning of the 



initial model. (The couplings Qc and Qs in ^] correspond to the couphngs for the osp{2\2) 
and su{2) currents respectively.) 

As before, in the second stage we restore other kinds of disorder by examining the quan- 
tum numbers of the remaining degrees of freedom. The residual currents are bilinears of 
fields in the [0, |]°*^ thus their quantum numbers follow from eq. (HTI) . Since A°g*^j = 
and the [8]°^^ has A = 0, for any A^ the dimension zero logarithmic field $[§] is the most 
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relevant operator. Therefore, the additional kinds of disorder should correspond to the 
logarithmic perturbation: 

f (fx 

S = Sospi2\2)^2N +9s J ^ ^[8] (58) 

where 5'osp(2|2)_2jv formally represents the conformal field theory with the current algebra 
symmetry. For N = 1 the detailed study of the effective action indeed leads to the same 
conclusion [itI. [s^]- 

For the special case of = 1, osp(2|2)_2 has a free field representation with the same 



content as for gl{l\l)k, so that Sosp{2\2)^2 has the free field form in eq. fl22l) [35l. |47|. The 



explicit form of the perturbation was given in |35||: 

$[8] = 4x'x' (cosh((0i - ^')/V2) + cosh((0^ + ^^)/V2)')+4{d,x%x'){x'x') cosh(v^0i) 

(59) 



54] 



For general the free field representation requires more fields 

Since the additional forms of disorder correspond to the above logarithmic perturbation, 
as we argued above, the exponents for the SQHT should be contained in the osp(2|2)_2Ar 
theory. For the 1-copy^ theory, this is more transparent using the gl{l\l)2 embedding since 
Tgi{i\i)2 = ^osp(2|2)_2[35]. The vector representation of osp{2\2) corresponds to the 5f/(l|l)2 
fields ^(1,0) and $(o,i) with scaling dimension 2A(i o) = 1/4, and this determines the density 
of states exponent p{E) ~ E^^"^ , since 1/7 = A/(l — A) with A = 1/8. The remaining low 
dimension fields are $(2,1) and $(1,2) with A = 5/8, —3/8 respectively. The other fields have 
dimensions which differ by an integer from the fields considered thus far. The A = 5/8 field 
determines the correlation length exponent for percolation: Upcrc = (2(1 — 5/8))^^ = 4/3. 
(In the osp(2|2)_2 description, the A = 5/8 field is a descendant of the field [±1, |]°*^ with 
A = —3/8.) Both of these exponents agree with the exact results in [^j- Note that the 
c = minimal model field with A = 1/3, which determines the localization length for 
self- avoiding walks, is not contained in the spectrum, which is consistent with observations 



made in 
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56[. Our proposal appears to be consistent with other observations made in 
56l |. since, because of the logarithmic perturbation, the critical point is not strictly speaking 
a conformal current algebra, even though it has some of the same exponents. A further check 
will be given in the next section based on the multi-fractal exponents. 

On the other hand, the above spectrum does not contain the full spectrum for the SQHT 
proposed in[56(]. The latter was based on a specific mapping of the lattice model to a 
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Coulomb gas. This is not necessarily inconsistent with our results, since the above spectrum 
is based on the simplest closed operator algebra and it contains the main exponents that 
are known to be physically meaningful. It seems likely that other twisted x sectors and/or 
bosonic cf) sectors could be consistently added to our theory, however we have not studied in 



detail the possibility of obtaining exactly the partition functions proposed in [56 1. 



VIII. MULTI-FRACTAL EXPONENTS. 



Generalities. 



In order to study multi-fractality in the density of states, we add an energy term in the 
action corresponding to H ^ H ~ S: 

S£ = J^^^ (^-^+ + ^-^+) (60) 
The density of states operator is then 

='^_^+ + ^_^+ (61) 

Multi-fractal properties refer to disorder averages of q-th moments of p, pi, and are simply 
related to wave-function ip multi-fractality since p = ip'^ip. Properly normalized quantities 
are 

{jdHp{x)Y ^ ' 

where here p represents (p) at fixed disorder. At a critical point P^'^^ scales as 

Pfa) L~^i (63) 

where L is the system size. The exponents Tq are related to the scaling dimensions of 
operators as follows: 

r, = r, - qT, + 2{q - 1) (64) 

where Tg is the scaling dimension of p"^ in the effective disorder averaged theory. 

For both the QHT and SQHT, there is a regime at low q where Tg is quadratic in q. Since 
Ti = and To = —2, in this parabolic regime Tg is characterized by a single parameter oq: 

Tq = {2- ao)q^ + aoq-2 (65) 
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Since the qTi term is simply a matter of normalization, it is meaningful to define 

f , = r, - qT^ = (ao - 2)g(l - q) (66) 
For the purpose of comparing with numerical simulations, one can perform the Legendre 



transformation 



55| 



dr 

/(a) = aq-Tg, " = i^'^) 



One finds q = {a — ao)/ (2(2 — ao)), which leads to 



^/ X (a — aoY , ^ 



-ao 



The parameter ao determines the typical density of states exp(log(p)) ~ L~ 

A common belief expressed in the literature is that in the supersymmetric method for 
disorder averaging, g-th moments can be calculated in the A^-copy theory for any N greater 
than q and one should obtain a result independent of A^. This can be proven by simple 
manipulations of the path integral, and this property is manifested for instance in the A^- 
independence of RG beta functions [3 7]. For the QHT (SQHT), it follows from properties of 
the gl{N\N) {osp{2N\2N)) BRST symmetries of the theory before and after disorder aver- 
aging. This "copy-symmetry" is analagous to replica-symmetry. In the scheme proposed in 
this paper, this copy-symmetry is broken in the first stage of the RG flow since gl{N\N) is 
broken to gl{l\l) for the QHT and similarly for the SQHT, as emphasized in section VII. 
Thus, the g-th moments will depend on the number of copies N one started with. Since 
this is contrary to expectations, at this stage it must be viewed as a working hypothesis. 
However we can provide some justification for this copy-symmetry breaking. First of all, the 
breaking of symmetries in the flow to a low energy fixed point is a common phenomenon, 
and in higher dimensions is more the rule rather than the exception. In the present con- 
text one must bear in mind the Mermin- Wagner theorem which prohibits the spontaneous 
breaking of symmetries. However the breaking of SU (2) SU (2) to SU{2) based on spin- 
charge separation in the Id Hubbard model, which is the prototype for our RG scheme, is 
well- under stood and known not to violate any theorems: the symmetry is broken in the RG 
flow but not spontaneously and there are no Goldstone bosons. Secondly, in spite of the 
powerful map to percolation for the SQHE, it remains unknown how to obtain the multi- 
fractal spectrum that has been found numerically, i.e. Fg ^ g(l — g)/8 for continuous q from 
this map. (See below.) Although it may just be a matter of time before this is eventually 
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understood, it could instead suggest that some basic assumptions are incorrect. The fact 
that we can obtain this result quite easily in the copy-symmetry broken fixed point is a 
positive indication. Thirdly, on the face of it, unbroken copy symmetry seems incompatible 
with the phenomenon of multi-fractality termination, whereas our scheme in fact relies on 
it to fix N, as we now describe. 

For q greater than some critical value Qc, Tq is no longer parabolic. This phenomenon of 

this paper, 
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59|. 



multi-fractality termination is thought to be distinct from the considerations o: 
i.e. it is a separate issue unrelated to the RG fiow of the disorder couplings 
For the QHT, 2 < gc < 3, whereas for the SQHT, 3 < gc < 4. For the SQHT the map 
to percolation was used to obtain only for the integer values q = 1,2,3 in [4^ and the 
argument breaks down for g > 3. 

If one wishes to compute then for all of q continuous and less than qc in a single, 
fixed theory, then this should be possible in the N copy theory with fixed to be the 
largest integer less than qc, i.e. = 2 for the QHT and = 3 for the SQHT. We wish to 
emphasize again that this new approach gives results that are different from the spectrum 



of multi-fractal exponents studied in 13|, l58|], and also differs from the calculation in [37|]. In 
the rest of this section we show that these assumptions appear to give results in very good 
agreement with numerical work. 



B. QHT 

The energy operator corresponds to the gl{l\l) fields 

p = $(i,o> + $(o,i> (69) 

We thus first need the gl{l\l) content of p'^ for integer q. Using 

{h,j) ® ((1, 0) © (0, 1)) = {h + © {h,j - 1) © {h,j + l)(B{h- (70) 

one sees that the above tensor product involves /i+j|ncw = h+j±l and /i— j|ncw = h—j±l. 
Therefore p'^ contains the representations {h,j) with —q < h + j < q, —q < h — j < q. 
Examining A^/jj) one finds that the most relevant operator in p'^ has h = 0,j = q. As 
explained above, we now set A^ = 2 = fc, and the field ^(o,g) has dimension Fg = 2A^o,g> = 
g(2 — g)/4, which gives 

f, = ^ (Tl) 
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i.e. ao = 9/4. We are implicitely analtyically continuing the 7" quantum number to continous 



4lj, ao = 2.26 ± .01, and in |42|, 



q. This agrees very favorably with the numerical results in 
ao = 2.260 ± .003. 

It needs to be emphasized that the copy-breaking feature of the RG flow in our scheme 
implies that Fi = 1/4 in the N = 2 copy theory is not supposed to equal the usual = 1 
result, i.e. Fi = 0. 



C. SQHT 

For the SQHT the density operator corresponds to the osp{2\2) field: 

p = (72) 

We need the quantum numbers of p'^ for q = 2,3. We have already considered g = 2 in eq. 
f HT]) . Since the [8]°'^^ can be viewed as [|, ® [~|) \]°^^, taking one more tensor product 
and using the rules in [l^ one obtains 

[0, |] ® [0, 1] ® [0, 1] = [0, |] © 3[0, i] © 2[i, 1] © 2[-i, 1] © [1, i] © [-1, 1] (73) 

where [b, s] refers to [b, s]°^^. Examining the conformal dimensions ^"^^g^, one finds that for 
q = 1,2,3 the most relevant operator in p'^ corresponds to [^y^, ^]'''^^- Fixing = 3, i.e. 
k = —6, as suggested above, one then has F^ = 2A°^^i ^ = g(2 ~ q)/8 which gives 

r, = ^ (74) 

In this case one has ao — 2 = 1/8. Again Fi = 1/8 in the A^ = 3 theory is not the same as 



the well 
In 



mown Fi = 1/4 of the A^ = 1 theory (see section VIIC). 
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431] it was found numerically that ao — 2 ^ 1/8 in the parabolic approximation 
(within a few percent). However it was also observed that in comparison to the QHT, for 
the SQHT there are more marked deviations from parabolicity. If one takes into account 
the non-parabolicity and simply defines ao from the maximum of /(a), then one obtains the 
result ao — 2 = 0.137±0.003. These observed deviations from a parabolic regime could have a 
number of explanations in our model. It could be due to the effects of logarithmic corrections 
due to the ^[s] perturbation. It could also be due to the large number of operators in eq. 
3]), where we took only the most relevant. 
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IX. LOCALIZATION LENGTH EXPONENT. 



In simulations of the network models, one needs to tune to a critical point by adjusting 
a parameter A, analogous to tuning to the critical probability Pc = 1/2 in 2D classical 
percolation. In our description, this should correspond to a term in the action 

SS, = I ^ Aa(x) (75) 

for some operator O^. If O^, has scaling dimension T^, then A has dimension 2 — Tu and the 
correlation length diverges as ~ (A — Xc)^'^ with u = 1/(2 — Tu). 

In contrast to the density of states exponents, we do not have arguments based on quan- 
tum numbers to identify the field O^. In the SQHT we know that T^, = 5/4. In the gl{l\l)2 
embedding in the 1-copy theory for the SQHT the field (9^ thus corresponds to ^*(2,i), which 



has the explicit form 35|] : 



$(2,1) = ^1/2/11/2 e^(^^^-^^)/^ + a?/,ai/, e^^Vv^ (76) 

Above, /ii/2 and o"i/2 are twist fields with conformal dimension —1/8 and 3/8 respectively. 

Lacking a first-principles identification of Oi, for the QHT, we can only give plausible 
values based on the spectrum of dimensions in our model. It was understood long ago 
that one must consider at least N = 2 copies, since the exponent describes criticality in 
the conductance. The latter is related to a product of retarded/advanced 2-point Green 
functions, and one needs separate copies for retarded verses advanced. We have already 
used the N = 2 copy exponents to explain the multi-fractality in the density of states. 
Therefore, it seems likely that one needs to consider N > 2 copies. Let us therefore double 
the number of copies one more time and consider the = 4 theory. As explained in 
section VII, when N > 2 there are potentially dangerous perturbations corresponding to 
the operators $(0,2) and $(i^„i) since ^(o)(4) is no longer the most relevant operator. The 
possibility of such additional perturbations can be further investigated, however since this 
is beyond the original scope of this paper, let us simply assume this issue is not important, 
and consider the theory 5f/(l|l)4 plus the logarithmic perturbation ^'(o)(4) as in eq. ( 15^ . The 
operator that most closely parallels Oi, for the SQHT is ^{n,n-i) with conformal dimension 

_ 2A(A-1) + 1 
A(7v,jv-i> - ^ (77 j 
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in the A^-copy theory, and leads to u = N'^/{2N — 1). At = 4, the operator $(4,3) has 
A = 25/32 and has the form 

ih _ ,, 77 ^i(4(/)i-3<^2)/2 I 2 —1 ^i(3</.l-2(/.2)/2 /yoX 

^(4,3> - /^l/4Ai3/4e ^ ^ ^" + cri/4'73/4e ^ ^ ^" (78j 

Here the twist fields /ii/4 and cr^^^ have A = —3/32 and 5/32 respectively. The above fields 
are local since one identifies the above operator with Oj,, then this gives 



the exponent u = 16/7, which is within 2% of the numerical results 2.35 ± .03 601] and 



2.33 ± .03 61[. (For a survery of the various methods see 221.) The value 16/7 ~ 2.29 is also 



consistent with the experimental measurement 2.3 ±0.1 

X. CONCLUSIONS 

In summary, by carrying out the RG fiow in two stages and using a new form of super 
spin-charge separation, we argued that the disordered Dirac fermion theories for the QHT 
and SQHT are described by logarithmic perturbations of the current algebras gl{l\l)j\f and 
osp{2\2)-2N where A^ is the number of copies. The explicit forms of the resulting actions 



were constructed using the recent results in [351]. We also argued that the logarithmic 
perturbations do not modify the exponents of the current algebra theories, however they 
generally lead to logarithmic corrections to correlation functions. 

The unconventional outcome of our RG scheme is that the A^-copy symmetries implicit in 
the gl{N\N) symmetry of the theory before disorder averaging are broken in the fiow to the 
low energy fixed point in the first stage since the SU {N)q is gapped out in the fiow. This leads 
to A^-dependence of the g-th moments, contrary to common expectations. Although this 
may seem problematic, we gave several arguments in favor of it in section VIII, and showed 
that it leads to a computation of the multi-fractal exponents in the parabolic regime, which 
agree very favorably with known numerical results. We also speculated on the localization 
length exponent for the QHT, and suggested one needs to consider A^ = 4 copies. After 
making some plausible assumptions, we were led to suggest the value u = 16/7, however 
more investigations of this proposal are clearly necessary. 

The QHT and SQHT are the smallest members of two of the main classes of disordered 
Dirac fermions 11], and the chiral GUE class was already solved in [16]. It would be inter- 



esting to investigate if the methods in this paper could be extended to the other classes as 
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well. 
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XII. APPENDIX A: SUPER-CURRENT ALGEBRAS. 

Let g denote a finite dimensional super Lie algebra and {J"} it generators. Each generator 
J" can be assigned a grade [a] = for bosonic generators and [a] = 1 for fermionic ones. 
The super Lie algebra can be presented as 

JV^ - j^j'^ = ffr (79) 

In the super current algebra at level k, denoted g^, the above generators are promoted to 
fields satisfying the OPE 

r{z)j\0) ~ ^r]''' + ^ffriO) (80) 

The currents have the mode expansion J"(z) = Xln'^n-^"""^ ^'^'^ ^^^^ modes Jq satisfy 
eq. (EH]). 

An important construction is the Sugawara stress tensor built on the casimir: 

T = kJ]JV^ (81) 

a 

where the coefficient k, is fixed by the requirement T{z)J°'{0) ~ J"(0)/z^. The conformal 
scaling dimension of a primary field in the representation r of g, is given by KC2 where C2 
is the quadratic casimir for r. 

Our conventions for the level k are based on the free field representations of these algebras 
in terms of the fields 'ip±,P±. Let = {ip±, P±) denote 2 component fields and 
a = 1,2, ..,N the iV-copy version. Let us arrange all these fields into 2N component fields 
a = 1,..,2N and let [a] be the grade. A complete basis of currents is defined in eq. 
(HM . OPE's can be readily computed from the OPE's of the in eq. (fT3l) and will serve 
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as our defining relations for tlie level k = 1 super current algebra. For instance, 

k 1 



where k = 1. There are similar relations of the form H{z)J±{Qi) ^ J±/ z and J+(2;) J_(0) ~ 
k/ z"^ + H/ z. We take the above OPE's as the defining relations of osp{2N\2N)k for general 
k. The H"'^ form a closed subalgebra which defines gl{N\N)k- 

Generally the Lie super-algebra g can be decomposed into its bosonic generators g*^°^ 
and fermionic generators g^^-*, g = g'-'^-' © g*-^-*. The bosonic generators form a closed sub- 
algebra, and the fermionic generators fall into representations of g'-^^ For osp{2N\2N), 
g(°) = so{2N) © sp{2N) and g*^^) corresponds to the representation (2A^, 2A^). The dimen- 
sion of osp{2N\2N) is thus 8A^^ and its rank is 2N. 

The algebra gl{N\N) is psl{N\N) ®u{l) ®u{l) where psl{N\N) is denoted A{N-1,N- 
1)/Z in 631]. The bosonic subalgebra is g'-^-* = sl{N) (Bsl{N) ©u(l) ©u(l). It has dimension 



4A^2 and it's rank is 2N. 

For the current algebras gl{l\l)k and osp{2\2)k we adopt a more specific notation. Again 
our conventions for the levels are based on the = 1 copy of the fields ip±, (3±- Define: 

S± = ±i^±f3^, Si = ^t/?t (83) 
These currents satisfy the OPE's 

H{z)H{0) ~ ^ 

2k 4 

z"^ z 

J{z)S±{0) ~ ±-S± 

z 

H{z)S^iO) ~ T^Si 
Mz)%{0) ~ (84) 

z 

z^ z 

-\ + -{H + J) 

Z'^ z 



J{z)J{Q) 


k 

r\_/ 

Z^ 


J{z)Mo) 


2 , 

~ ±- J±, 

z 




z 


H{z)S±iO) 


~ ±-s±, 

z 




2? 

z 

S^{z)S40) 
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Since sp{2) = su{2), the bosonic subalgebra generated by H, J, J± corresponds to su{2)Q)u{l) 
where the J, J± are su{2) generators. Rescahng J 2J, J± ±.2^/2J± one finds they 
satisfy the su{2) current algebra at level —k/2. Thus, osp{2\2)k has an su{2)_k/2 subalgebra. 

Our conventions for gl{l\l) are taken from the subalgebra of osp{2\2)k generated by 
S±,H, J and are presented in eq. (|T71) . The stress tensors and scaling dimensions of primary 
fields are given in the main body of the paper. 

Note that letting H ^ J, J ^ H, S± ^ ^S±, the new currents satisfy fc. Thus 

there exists an automorphism of gl{l\l)k which fiips the sign of k, and was used in 35|. 
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